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^_^ Abstract 

■•^ . In an extreme binary black hole system, an orbit will increase its angle 

I of inclination (t) as it evolves in Kerr spacetime. We focus our attention 

?-H on the behaviour of the Carter constant (Q) for near-polar orbits; and 

, ^■*J develop an analysis that is independent of and complements radiation 

reaction models. For a Schwarzschild black hole, the polar orbits represent 

^H the abutment between the prograde and retrograde orbits at which Q is 

^ at its maximum value for given values of latus rectum (Z) and eccentricity 

^\ (e). The introduction of spin {S = IJj/M'^) to the massive black hole 

OO causes this boundary, or abutment, to be moved towards greater orbital 

^~~^ inclination; thus it no longer cleanly separates prograde and retrograde 

\l^ orbits. 

[~ — . To characterise the abutment of a Kerr black hole (KBH), we first 

f^ investigated the last stable orbit (LSO) of a test-particle about a KBH, 

^^ and then extended this work to general orbits. To develop a better un- 

^~~^ derstanding of the evolution of Q we developed analytical formulae for 

^ Q in terms of I, e, and S to describe elliptical orbits at the abutment, 

•^ polar orbits, and last stable orbits (LSO). By knowing the analytical form 

/\ of dQ/dl at the abutment, we were able to test a 2PN flux equation 

^ for Q. We also used these formulae to numerically calculate the di/dl 

of hypothetical circular orbits that evolve along the abutment. From 
these values we have determined that di/dl — — { 122. 7S — 365"^ ) P^^" — 

63/2 S + 35/4 S^) ^^/^ - 15/2 Sr'^''^ - 9/2 Sr'^''^. By taking the limit 
of this equation for / — > oo, and comparing it with the published result for 
the weak-fleld radiation-reaction, we found the upper limit on 9t/SZ for 

the full range of I up to the LSO. Although we know the value of dQ/dl 
at the abutment, we find that the second and higher derivatives of Q with 
respect to I exert an influence on di/dl. Thus the abutment becomes an 
important analytical and numerical laboratory for studying the evolution 



of Q and t in Kerr spacetime and for testing current and future radiation 
back-reaction models for near-polar retrograde orbits. 

1 Introduction 

In his landmark work of 1968, Brandon Carter derived a new constant of mo- 
tion that pertained to orbital motion in the gravitational field of a Kerr black 
hole (KBH) [I]. In due course, this constant became known as the Carter con- 
stant, which joins the set of important constants of motion: orbital angular 
momentum {L^, z-axis projection), orbital energy (E), and finally the Carter 
constant (Q). These constants of motion can be developed rigorously from the 
Hamilton- Jacobi equation |T1[31IS|. 

An extreme mass ratio binary black hole system is composed of a secondary 
object (which may be a compact object (CO) of several solar masses) in orbit 
around a primary object (which is a massive black hole (MBH) of several mil- 
lion solar masses). Extreme mass ratio inspirals (EMRIs) are expected to emit 
gravitational wave radiation (GW) of sufficiently high energy and in the appro- 
priate frequency band for detection by the Laser Interferometer Space Antenna 
(LISA) to be feasible jU [SJ IHl [Z] • The emission of GW causes the constants of 
motion to evolve, which in turn affects the GW power spectrum. Therefore some 
useful methods have been developed to describe this evolution. For example, 
the quadrupole formalism |H1 El [M E] and the Teukolsky equation [HI ESI E] 
have yielded important results. The analytical description of the evolution of Q 
has been more difficult to achieve than it has for the other two constants of mo- 
tion [TS] ; although the use of the Teukolsky equation has shown great promise 
[ini ESI H |51 17] in this endeavour. 

As the CO inspirals, the gravitational radiation reaction causes the value of 
Q to change |inilI71EElESllllEIll2nilIIllIll2aiI|. Therefore a non-equatorial 
orbit lists as its angle of inclination, l, increases with respect to time; a near- 
polar prograde orbit becomes polar, and ultimately retrograde [53J EZ|- Such 
listing behaviour of an inclined orbit has been studied and confirmed using the 
most current Teukolsky-based fluxes |52] . It is our goal to develop an analytical 
and numerical methodology for testing and improving radiation reaction models 
for predicting orbit listing and inspirals for near-polar orbits. 

Our interest lies in studying KBH systems; yet, the Schwarzschild black 
hole (SBH) is an important datum. An infinitesimal amount of spin angular 
momentum {SS <C 1) may be imparted to an MBH such that, for practical 
purposes, it can be regarded as an SBH (by virtue of its minuscule effect on 
the surrounding spacetime); and yet, the spherical symmetry of the system has 
been broken and a z-axis defined. Then an SBH can be considered to have 
a prograde or retrograde inclined orbit; and the set of polar orbits define the 
abutment, at which Q will be at its maximum value (Q is non-negative for any 
bound orbit), for given values of /, e, and S (ceteris paribus). 

In section [2] the motion of a test-particle in an inclined orbit is analysed 
from first principles |24[ |25| |26| |27] to yield the effective radial potential and 



an analytical expression of L^ for a last stable orbit (LSO) . In section bl we 
continue our analysis to find the roots of the effective radial and polar-angle 
potentials and use them to derive analytical expressions for E and X"^ (where 
X = L^ — SE). The concept of the abutment is then refined. In section Ul we 
derive a set of critical formulae that express Q at the LSO, along the abutment, 
and for the set of polar orbits. The interrelationships between these formulae axe 
examined and a map of admissible values of Q, with respect to / and e is drawn. 
In section [5J this map is used to better understand the path in the Q-l plane 
that is followed by an evolving circular orbit. We demonstrate the importance 
of the first and second derivatives of Q (on the abutment) with respect to I for 
understanding the rate of change of t as the orbit lists. 

We shall use geometrical units by setting the speed of light and gravitational 
constant to unity (i.e. c — 1 and G — 1); therefore, mass-energy is in units of 
time (seconds). In addition, many of the parameters we use will be normalised 
with respect to the mass of the black hole (M) or with respect to the test- 
particle mass (m) . In appendix A, the symbols used in this paper are tabulated. 
By emphasising normalised variables, the analytical equations and numerical 
formalism are much better handled. 

2 The Motion of a Test-Particle in an Inclined 
Orbit 

A sound mathematical analysis can be made on the assumption that the sec- 
ondary body is of infinitesimal mass (i.e. a test-particle). In such a case, the 
background metric of the MBH dominates. In the case of EMRIs, the small 
ratio of the CO mass to the MBH mass (77 ^ 10^^) warrants our use of idealised 
test-particle calculations [U [2SJ |5J [7] . 

2.1 Basic orbital equations 

We begin by considering a test-particle in orbit about a KBH of arbitrary spin, 
S, for which the four-momentum can be given the general definition |29| . 



P. 



-mE, m-^ 

A 



dR 



, mMLg, mMLz 



where 



and 



A = [R^ -2R 
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(1) 
(2) 



S = 



R' 



S^cos^ [9) 



(3) 



Unlike the analysis in Komorowski et al. ^^, we shall use normalised variables 
at the outset and oS^er a more thorough treatment. Because we are now con- 
sidering inclined elliptical orbits, one cannot simplify the four-momentum by 



setting Lq = 0. But by knowing the Carter constant in terms of normalised 
variables (i.e. obtained by dividing through by mM), 

Q - ^°'' i^jf' +11 + cos^ {&) S' (l - E^) , (4) 



sin^ (9) 

one can obtain the component of orbital angular momentum, L, projected upon 
the equatorial plane of the KBH, 



cos2 (9) LI 



(9)~S'^{\-E-^), (5) 



V sin2 [9) 

and substitute it into the expression for the four-momentum: 
Pi = (6) 

' mE.mt {§) ,mM^Q - ^^^^ - cos^ (9) S^ (l - E^),rnML, 

The invariant quantity, 

PP ^ P,Psg'^\^^^,^ = -m\ (7) 

is calculated tensorially using the inverse Kerr metric (see Appendix B) and 
used to develop the radial orbital equation for a test-particle: 

2(dR\ _ ( £,2\ n4 , o n3 



S^l^l = -{l-E']R'' + 2R 

'lI + s^[i-e^ )+q)r 

2 
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2 Ul,-Se'^ +q]r-QS^ (8) 



By setting Q — 0, equation (Isl reduces to the equation for an equatorial orbit 
(see EO]) 

At the radial turning points, dR/df = 0. Equation (Isl becomes: 

= R^~2 ^ 



1-^2 
L2+S'2fl-^2 )+Q)i?2 



1-E^ 

l^-seV + q)r ~2 

2^ 4 ^ + ^^. (9) 



In the limit 5 — ?> (set S — while retaining a non-zero value for Q) 
equation (|9l becomes 

^3 (lI + q)r' (lI + q)r 

Thus, the square of the total orbital angular momentum, L'^ ^ Ll + Q, confirms 
that, for the specific case of an SBH, Q represents the square of the component 
of angular momentum projected on x — y plane of the coordinate system (see 
Appendix B in Schmidt |30| and Appendix C in this paper for a more detailed 
treatment) . 

Some important research [211 I31j has been performed by working with the 
orbital inclination angle, t, instead of Q; but in our study, the value of Q will be 
taken as a system parameter. If Q is chosen to be zero, then the orbital plane 
coincides with the equatorial plane of the KBH (i.e. l — and 6 = 7r/2) and for 
a test-particle in a polar orbit (i.e. i = f and < < tt) L^ must vanish. The 
choice of working directly with the Carter constant, Q, as a system parameter is 
consistent with the approach taken by Carter jT] and more recently emphasised 
by others [31 [3 [52 . 

2.2 Effective Radial Potential 

To proceed, we use a version of equation (|8|, which is quadratic in, E, 

-r(r^ + RS^ +2sAe^ + A RL^SE 
+R{R~2)(q + lI)+R'^(r'^ -2R + S^)+QS^ ^ 0(11) 



The roots of this equation can be used to determine the effective potential of 
the test-particle (V±): 



v± = (2rlJ±\/rza] 



-1 



X [R[R^ + S''R + 2S^]) (12) 



where 



Z = R''+{Q + S-' + Li]R-' + 2R'S^ + S^RQ + 2QS' 



When the last stable orbit (LSO) is reached, E corresponds to a local maximum 
of V+ closest to the event horizon. Therefore one calculates the derivative of V+ 
with respect to R and equates it to zero, i.e. 

dV+ ( ^ 



± =- 2RSL^ fa R^ + S'^\ \/rzK + R^Z^L 



dR \ 



R^ + S^R + 2S^]Z; 



I rVrZA(r^ + S^R + 2 sA 



(13) 



where 



and 



Zi ^R^~3R'^ + S^R^ ~3R'^S^ + 6S^R-2S'^, 

Z2 = -R'^ + R'^Q- (2S^ + 3q')r^+(2QS^+AsAr^ 



2.3 Orbital Angular Momentum at the Last Stable Orbit 

We can now develop an equation for the L^ of a test-particle in an inclined orbit 
about a KBH (in a manner similar to that described in [29 ) and extend the 
concept to general orbits. It should be noted, the value of L^ considered here 
is not valid for general orbits, but pertains to the LSO. After eliminating the 
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square root in equation (13 1 to yield a new equation that is quadratic in, L 
the solution is found to be: 

Ll = {-i?8 + (3 + Q) R' + (-2 52 - 6 0) i?*^ 
(-6 + 2 Q) 5' + 9 g) i?5 + (-54 + (-10 g + 12) 52) i?4 

+ ((-5 + Q)S^ + Q QSA R^-6 S'^R^Q 

+5 5^i?g-2S'^g 

±25(3 R"^ + sAaMr^ - R^Q + 3 R^Q + Q^sA } 

N -1 

R^(r^+9R-4:S'^-6rA] ,(14) 



which provides a relationship between L^ and R (which here, represents the 
pericentre radius) of the test-particle LSO. This result is independent of whether 
one begins the calculation with V+ or V- . 

One can now plot Lz with respect to the value of the pericentre (the point 
of closest approach, Rp) for an LSO for the cases where S — 0.0 (Figure ll|, 
S = 0.5 (Figure [2]) and S = 0.99 (Figure [s]). The values of L^ calculated for an 
SBH are plotted in Figure [l] for the range of Q values 0.0 to 16.0. 



For an SBH, the prograde (plus) and retrograde (minus) formulae for X^ 
(equation ( [l4| ) are reflections of one another about the R axis (Figure [11. But 
when iS* > 0, the plus equations are pulled below the R axis and this symmetry 
is lost (Figures [2] and Isl. There now exists a set of retrograde LSOs which 



are governed by the plus form of equation (14). The importance of this fact is 



revealed when we set the quantity beneath the square root in equation ( 14 ) to 



R" -QR'^ + SQR-' + Q^S' = 0. (15) 

The polynomial describes the boundary at which the plus and minus equations 
for L1 are equal and it offers an insight into the behaviour of Q for LSOs that 
are nearly polar (t ~ 7r/2). If 5 == 0, then for an ehiptical LSO [33l|TTl|^. 

Rp = 2{3 + e)/{l + e). (16) 



By substituting equation ( 16 I into equation ( 15 1 and solving for Q, one obtains 



Q = 4(3 + e)'[(l + e)(3-e)]"\ (17) 

This result applies to LSOs at the boundary and specifies an upper limit on Q 
for orbits around an SBH. Now we must develop these ideas for general orbits 
about a KBH that have not yet reached their LSO. 



3 Analysis of the Trajectory Equations 

3.1 Introduction 

There exist four trajectory equations [2 [Ml [3U] in two categories: 

category (a) 

(those that are periodic in radius, R, or polar angle, 9) 

t^ = ±^VniR) (18) 



S^ = ±V/^ (19) 

category (b) 

(those that are monotonically increasing in azimuthal angle, <p, or coordinate 
time, t) 

t% . V,^ (21) 



See Appendix B to see equations for functions Vr {R), Ve, V^, and Vt- We have 
already developed equation (18) in section 2.1 (equation (Isl). And equation 



(19) can also be developed by a similar method (see Appendix C). 
One obtains, viz. equations (18) and (19), the following condition 



dR 



d6 



Vr (R) 



Veie) 



(22) 



on the geodesic of the test-particle, which is a general form of the equation 
specified by Schmidt (equation (16) in [30). Given equation (19) [34J, one can 
find the proper time of the orbit: 



T^de 



Vg 

"-' R^dR 



'Vr 



S' 



' cos2 (6*) 



d0. 



(23) 



where the integral has been split into its separate R and 9 integral terms viz. 



equation (22). The same result is found when starting with equation (18) in- 
stead. Two other important integrals that can be calculated for coordinate time 
and azimuthal angle are given in Schmidt (equations (14) and (15) [3U]) in which 
a detailed analysis is made on the basis of the Hamiltonian. Equation (23 1 can 
be solved to yield elliptic functions [.2DJ; therefore, the roots of Vr and Vg con- 
tain information necessary for deriving analytical formulae for Q in terms of I 
and e (for given S) and t as a function of Q. This will be shown in sections 
andlD 



3.3 



3.2 Roots of the Radial Equation 

We introduce X ~ Lz — SE and convert equation ([9| to the form: 

= - - «' 





'{^--') 


(A- 


f 5*2 + 2SEX + q) i?2 




(i-fi^) 


2 {X 

(' 


' + Q)R QS^ 
-£2) (l--^') 



(24) 



This substitution is consistent with the approach in [29; and that undertaken by 
Glampedakis and Kennefick [20]; and it will help us to derive the latus rectum 
of the LSO. Analytically, the use of X^ in this case offers an advantage over the 
use of L^ . 



3.2.1 Elliptical Orbits. 



By finding the four roots of equation ( 24 1 one can derive analytical formulae 



for X and E, in terms of e, I, Q, and S, which apply to general orbits (and 
are not limited to the LSO). Although the roots are easily obtained in terms of 
the constants of motion: L^, E, and Q; they are complicated and as such not 
helpful. To simplify the analysis, we assume a priori that an inclined orbit can 
be characterised by an eccentricity, e [3S]. Therefore the radius of the orbit at 
its pericentre is described by 



/ 



^ 1 + e' 

and correspondingly, at its apocentre 

I 



1-e 



(25) 



(26) 



To proceed, consider an expansion of the four possible roots {r^ < r^ < rp < Va}'- 



- {rs + r4 + ra+ rp) R^ 

+ (^^3 + rsra + r^rp + r^Va + r^Vp + rpra) R^ 

- {riT^ra + r^r^^rp + r^rpra + r^rpra) R 

+r4r3rpra = (27) 



By equating the coefficients of the two polynomials in equations ( 24 ) and ( |27| 
one obtains the two independent equations: 



and 



I [2 r^r^ + r^l + rj\ X^ + Q 



r^r^f 
l-e2 



1-^2' 



(28) 



(29) 



which have been simplified viz. equations ( 25 1 and ( 26 ) . Let us solve equations 



(28) and (29 1 to obtain 



l~e' 



r3 



p (i - eA 



{q{i 



l-S^] +X^l] ±Jz. 



(30) 



and 



where 



r4 



QS^ (1 - e^) 
rsP (l - E^ 



(31) 






^4r2 



l-e2 



-2Xn + 2X^S^ \Q + X*l 



If Q = then selecting the minus sign in equation ( 30 1 yields r^ — 0; therefore 



the plus sign is the one taken as physically meaningful. For Q = 0, equation 
( 30 ) reduces to 



X^ (1 - e2 
P 1-E^ 



(32) 



which applies to an equatorial orbit. The value of r^ equals zero when Q = 



as can be seen in equation (31 1 



3.2.2 Parabolic Orbits. 

Parabolic orbits have importance to the empirical study of the interaction of 
stars with massive black holes (MBHs) [36 , 37 . For parabolic orbits both e = 1 
and E = 1. We refer back to equation ( 24 ) ; and set E = 1: 



= R-" 



1(^2 , 52 



{X^ + Q)R^-QS' 



2SX + Q]I{- 
1 



(33) 



There are now three possible finite roots {r^ < r^ < Vp}, which can be used in 
a new general equation (r^ is infinite in the case of a parabolic orbit): 



R^ - (^3 + ri + rp) i?2 
+ {riTs + r^rp + r^rp) R - rir^rp = 0. 

The pericentre simplifies (viz. e = 1) to become, 



(34) 



2' 



and correspondingly, at its apocentre. 



(35) 
(36) 
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Q{l-S^]+X^l]±^/Z. 



We obtain the solutions for the additional roots 

1 

r-^ — ^^ 

P 

and 



TA 



rzl 



where 



Zi = [l-S^\ Q^ -l{Sn^ -2X^l + 2X^S 



32 r2 



^2 o2 



3.3 Roots of the Polar-Angle Equation 



Let us focus on the denominator of the second term in equation (23 1, i.e 



XH 



(37) 

(38) 

472 



to derive an analytical relationship for the limits of integration, 9i and 6*2, from 
which one may determine t. We shall work with Ve in terms of L^, i.e. 



S' 



sin (6i)cos2 {d)de 



Q sin^ {6) - sin^ {0) cos2 {6) S^ (l - eA - cos^ {6) hi 



(39) 



Equation (23 1 is an elliptic integral; thus the limits of integration correspond 

(40) 



to the zeros of the denominator. By making the substitution 

u — cos {Q) , 
the integral, I, becomes 



I 



S' 



= S' 



-u du 



52 1 - £;2 m4 _ g + 2. 



^ 2 



52 1 - £2 



«i S^/l - &^{v? - /3+) (m2 - /3^y 



Q 



(41) 



for which the roots, /3± — cos2 {9±), can be calculated and used to determine 
{viz. cos (9) — cos (f — i) = sin (l)) the exact angle of inclination of an orbit for 
which the values of S, Q, and L^ are known, when working in Boyer-Lindquist 
coordinates. We will calculate l using 



sin^ (i) = - 



Q + Ll 



52 1 - ^2 _ 



+ L2 + 52(l-^2 -4Q52(l-£;2 



5*2 (1 - eA 



.(42) 



which differs from the approximation in |211 [5]. N.B.: in dealing with results 
that are first-order and third-order in S, one may consider the approximation 
to i to be reasonably close to equation (|42l |38| . 
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3.4 Orbital Energy and an Analytical Expression for X^ 

As outlined in [20] , the next step will be to develop a formula for orbital energy, 
E, in terms of e, I, and X'^. By referring to equations (24) and (27 1, this 
derivation proceeds by solving 



to yield 



Ti + rs+rp + ra = 2 [l-E' 



(43) 



E = ±- 



1-e' 



P 



Q[l-SA +XH + P(e^ + l-l) (1- 



for which we use the positive case. For Q = 0, equation (44 1 simplifies to 



E = 



•l- 



I 



X2(l 



P 



(45) 



which is the expression for E in an equatorial orbit. 

It is interesting to observe tha t for an inclined orbit around an SBH (5 = 0), 
the equation for E (equation (44l) reduces to 



E = 



\ 



(f-2(l + e)) (l-2{l-e) 



(46) 



which is the expression for orbital energy of a test-particle in orbit around an 
SBH, (Cutler, Kennefick, and Poisson (see equation (2.5) in [33j)). Further, this 
equation for E shows no dependence on Q. This property is expected since the 
orbital energy must be independent of orientation in the spherically symmetric 
coordinate system of an SBH. 



In the general case, we observe that E' is a function of X^ (see equation (44 1); 
and thus it is not in explicit form because X — Lz^SE. Therefore the use of X^ 

in place of ( Z/2 — SE j simplifies the analysis by avoiding an unending recursive 

substitution of E into the equation. Although one may derive a formula for E 



in explicit form, it is better to perform the analysis using equation ( 44 1 

To calculate an analytical expression for X^, we substitute equation (|44l 



into our original quartic (equation (|9|)) and evaluate it at either Tp or r^ (the 
two simplest choices of the four roots) to yield: 



(1 + e)'^ (p (52 - /) + 2 QS^ (1 + e^) + I (X^ + Q) ( 



l-i-e' 



(47) 



+2SXJl{X^ + Q){l-e-^f -QS-^{l-e^f + l^[l + e-^ -l)\ = 0. 
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By eliminating the square root, and solving for X^ in the resulting quadratic, 
one obtains: 



r2 ^5 + ZqQ ± 25* \J Z^Z^Zi. 



where 



lil{3-l + e^j -4S^l-e^f] 



Zr, = -2 1 



and 



)' 54 + 2f (2e4 + (2 - Pj e^+4- Pj S^ 
-P{3-l + e^y , 

Zy = S^l + ef + 1 (l - 2(1 + e)y 

Zs = S^l-ef + l{l-2{l-e)y 

Zg = (P + S^Q^ (1 -e^f + QP {3-1 + e^)) . (49) 



Xj. has a minus and a plus solution, which we will carefully describe in the 
next section. We have avoided the analytical difficulties that would arise by 
working with L^ directly. Indeed, the advantage of using X = L^ — SE is more 
than a simple change of variables, but rather an essential step in solving these 
equations. 

3.5 Prograde and Retrograde Descriptions of X^ 



The expression for X"^ (see equation (48 )) contains the square root, ±25* \/ ZjZ^Zq 



for which, Q is found only in Zg as a quadratic. Therefore it is easy to derive an 
expression for Q for which Zg — Q and thus determine where the minus and plus 
equations for X"^ meet or abut {viz. Zq = 0). This information is important for 



determining which form of equation (48) to use. We will call the set of general 



orbits for which Zq = 0, the abutment, to avoid confusing it with the result for 
the boundary between the plus and minus forms of L^ at the LSO. 

A prograde orbit has an Lz > 0. Correspondingly, when Lz < the orbit 
is retrograde; and if Lz = 0, the orbit is polar. When using the minus and 
plus forms of the equation for X^. (equation (WSJ) ) , one must recognise that X"^ 
governs all of the prograde orbits, the polar orbits, and the near-polar retrograde 
orbits up to the abutment; and Xl_ governs the remaining retrograde orbits. If 
one considers an SBH system then the abutment will be comprised only of polar 
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orbits (and it is only then that X'j_ cleanly separates the prograde and retrograde 
orbits). If 5' > 0, the abutment will always consist of retrograde orbits. 

For orbits governed by X^, X+ — —JX'^; but for those governed by X^, 
the choice of sign depends on the value of Q. The plot of X^ with respect to 
Q for a circular orbit with / = 6.25 about a KBH of S^ = 0.5 (Figure U| shows 
that for dX^jdQ to remain continuous over the range of real values of Q, the 
minus sign must be chosen when evaluating X_ for Q > Q switch- An analytical 
formula for QswUch can be found by solving X^ = for Q. The general solution 
is 

Qsn,^tch = P {l- ^') {I (l~ e2 - 3) + 2 (1 + e^) ^2^ "'(^O) 



and for an SBH 



^switch = /^ ( / - e^ - 3 ) . (51) 



For large orbits {I — > cx)) equation ( 50 ) can be converted to its asymptotic form 



by first factoring out I from each term to obtain, 

for which the denominator may be brought up to the numerator to yield 





In the limit as / — ^ cxd, 



C switch 



I 



= l + 3 + e'^-S^. (54) 



Equation ( 54 1 describes the locus of points at which L^ = SE (which is effec- 
tively constant for large I); therefore, Q switch does not describe a trajectory. 
We have developed two formulae: one for the L^ at the LSO (equation ( 14)) 



and the other for the X'^ of general circular and elliptical orbits (equation (48 1). 
For each, there is an expression that describes where the plus and minus forms 
are equal. For Z^, the boundary between the plus and minus forms is described 

by 

Rl~QR^ + 3QRl + Q^S^ ^ 0, (55) 

where Rp represents the pericentre. And for X^., the abutment is described by 

l'^ + S'^Q^{l~e'^Y + Qp(3-l + e'^) = 0. (56) 



As equations ( 55 1 and ( |56[ ) describe the boundaries (where the plus and minus 
forms are equal) that pertain to different quantities (L^ and X|.) they will not 
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in general coincide, 
obtains: 



If one substitutes Rp = Z/ (1 + e) into equation (55) one 



l^ - Q(l + e)/* + 3Q(1 + eyi^ + Q^S^l 



0, 



(57) 



which equals equation (56 1 when e = 0, (i.e. for a circular orbit). If S" = 
(SBH case), then X^. ~ L\ and the two boundaries must be identical. If 
one substitutes Z = 6 + 2e into equation (56 1 and solves for Q, then the same 



expression as in equation (17 1 is obtained 



4 The Characteristics of the Carter Constant Equa- 
tions and the Domain of the Orbital Parame- 
ters 

4.1 Introduction 

In describing an arbitrary orbit, it must be recognised that each parameter (I, 
e, and Q) has a domain. The value of e lies between 0, for a circular orbit, and 
1 for a parabolic orbit. Although I has no upper limit, its minimum value is 
^LSO'i while Q, which is non-negative, does have an upper limit that depends 
on the size of I. The complicated interrelationships between these parameters 
can be better understood if we derive a set of analytical formulae to describe 
the behaviour of Q with respect to I and e. In the sections that follow, we shall 
examine the LSO, abutment, and polar orbits. A representative plot of these 
curves is shown in Figure [s] for e = {0.0, 0.25, 0.50, 0.75, 1.0} and a KBH spin 
of ^ = 0.99. 



4.2 Last Stable Orbit 

In |29| a new analytical formula for the latus rectum of an elliptical equatorial 
LSO was developed. We can perform a similar treatment for inclined LSOs; 
but the polynomial that results is of ninth order in I, currently making the 
derivation of an analytical solution infeasible. The use of the companion matrix 
[391 HUl [25] simplifies the numerical calculation of the prograde and retrograde 

hso- 

We refer back to equations (25 1, (26 1, (30), and (31); but because E will 



equal the maximum value of V+ (closest to the event horizon) we can specify 
rs = Tp as an additional condition |26| , I33| , [29]. Therefore the remaining root, 
r4, can be easily solved viz. 



Til^ 



QS' 



(l + e)(l-e2) 1-^2' 



(58) 



15 



to yield, 



ta 



QS^ (1 + e) (1 - e2) 



l-E^]P 



(59) 



Substituting the four roots (two of which are equal) into equation ( 43 ) yields 



QS^{l + e){l~^ I I 



1-^2/3 



1 + e 



^.2(l-E')-\ (60) 



Substituting the expression for E f equation (44l) into equation (60 I and cross 
multiplying, we obtain: 



l^ {1 + ey {I - ey ZioZu = 



(61) 



where 



and 



-(- 



Zw = [-31 - 2le + e^l] Xl + iQS^ + AQS^e - 2Qle + Qe^l - 3QI + f 



Zii = (1 - e^) (Ixl + Ql- QS^^ - P 



By setting Ziq = and substituting the formula for X^ {viz. equation (48)) 



one obtains a polynomial, p(l), in terms of / (to order 9), e, Q (to second 



order), and S (see equation (1031 in Appendix B). The companion matrix of 



p [l] provides a powerful method to calculate the values of I for both a prograde 
and retrograde LSO by numerically evaluating the eigenvalues of the matrix. 
Optimised techniques for solving for eigenvalues are available |41l H^ . Such a 
numerical analysis was performed to obtain representative values of Ilso- The 
corresponding values of Ilso 7 which we derived from Zn = (ceteris paribus) 
were smaller and thus not physically reachable by a test-particle. This result 
demonstrates that Ziq = is the appropriate solution. 

Because p (l) is a quadratic in terms of Q, an alternative way to analyse the 
behaviour of orbits as they approach their LSO is available. One solves for Q, 
analytically, to obtain: 



where 



Qlso — -7 [Zi2 — S Zi3^/Zi4J is Z15J , 



^12 



(62) 



[4 (e + 1) (l^ - (2e2 + e + 3) 1 + 2 {e + 1) {2e'' - e + 3)] S^ 
-P (e + if (2e^ + e^ + I {3 - e) + l) S•^ 
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Zi3 = (3-e)((e + l)'52-f(2e + 2-f)), 



and 



Z 



15 



-Up - (e + 1) (e^ - 2e + 3) [+ (2e2 - 4e + 3) (e + if)) 



We now have, through equation (62 1, the means to plot the value of Qlso 
with respect to Ilso- For a KBH, one may use Ilso in the domain Ilso, prograde-, hso, retrograd 
as the independent variable. 

4.3 The Abutment 



The roots of equation ( 49 1 allow us to calculate the value of Q along the abut- 



ment {Qx) in terms of the orbital values of /, e, and the KBH spin, S; and we 
obtain, taking the minus sign of the quadratic solution: 



Q 



X 



2S'2(l-e2)' 



('( 



/ - e' - 3 



Pil-e^-i] -U{l-e^fS'^ . 



(63) 



Equation (63 1 appears to have poles at S" = and e = 1; but we can demonstrate 



that it reduces to a well behaved function at these values of S and e. We first 
factor out the term, I (l — e^ — Sj , to obtain 



Q 



X 






( 



4(l-e2)2^2 



Iv \ '^ij~^'-^) 



(64) 



/ 



Equation ( 64 ) can be simplified by making a binomial expansion of its square 



root to yield, 

[3 (f- 6^-3 



Q 



X 



2(l-e2)'52 



^n-i + ^ '- — 

^■2 I r/r o „\ 



/([-e2-3y 



(65) 



{l-e'Vs^ 



I l-e^-3 






10 



:i-e2)'^2 



/ / - e2 - 3 
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Equation (65), simplifies to a power series in terms of (l — e^) S*^; therefore, 
only the first term will remain when S" = or e = 1. We now have a greatly sim- 
plified expression, which applies to elliptic orbits around an SBH and parabolic 
orbits about a KBH, i.e. 



Q 



X 



P {l-e^-3 



I' I 



Further, as ^ — > oo in elliptical orbits. 



4 1 



2 A9 



{1-eYS 



(66) 



(67) 



therefore, a similar binomial treatment may be performed on equation (64) to 
yield Qx = ^ + e^ + 3 for large /. This result is consistent with the fact that 
as I — )■ oo, the spacetime looks more Schwarzschild in nature. As in the case of 



equation (17), Qx (equation (63)) also defines an upper limit on Q for specific 
values of I and e. The points above the Qx curve are inaccessible; and this 
property can be shown by direct calculation. 

This result confirms the choice of the minus sign in solving the quadratic 



equation (56) for Qx- Recall that equation (48) applies to bound orbits in 
general and is not restricted to LSOs (as is the case for L^ (see equation (14))); 
therefore, the value of Qx can be evaluated for all values of / > lLSO,prograde', 
and it will apply to both sets of orbits governed by X^ or X"^ . 



4.4 Polar Orbit 

We shall consider polar orbits in this paper; although orbits of arbitrary inclina- 
tion are also important. The polar orbit (governed by X'^) is precisely defined 
by setting L^ = from which we obtain X'^ — S^E^; therefore, we can derive 
an analytical formula for Q of a polar orbit {Qpoiar) of arbitrary I and e. 

We have the formula for X^ in terms of I, e, and Q (see equation (48 1); 
and E can also be expressed in these terms (see equation (44 1). Therefore 
Xi — S'^E'^ = can be simplified and factored to yield: 



where 



A[BiQ-fB2] [CiQ-fC2 







(68) 



A = [3_(2e2 + 6)P + (e2 + 3)'/'-452(l-e2)' 



Bi = l^bi = 



f^jl - (3 + e^) f-i + 2 5*2 (1 + e^) t"^ - 2 5^ (l - e^fl 



-SUl-e^fr^ + SUl-e 



'fl-'}. 
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B2 = 1% = P |l + 2 5^ (1 + e^) l- 



-4S'2(l-e2)r 



SUl- 



'yr'Y 



Ci = rci 



l^il - (3 + e^) [-1 + 2 5*2 (1 + e^) r^ + 2S^{l~ e^f l 



and 



C2 



l^C2 



p{l- 

-4 5^ 



(l_e2)^~-4 + 54(^_g2)3^--_5|^ 



AsH-^ 



2 5*^ (3 
1 - e^) r-3 + 5* ( 



e2)r 



l-e-ft^]. 



The factor, ^, offers no physically meaningful results. It does not provide a 
solution for Q; and for ^ S" < 1.0 and ^ e ^ 1.0, we find that 3 < f < 4, 
which lies beyond the LSO. The factor, BiQ ~ PB2 = 0, yields the result 

Qpolar = I B2B^ . (69) 

And the factor, CiQ — PC2 — 0, yields the result 



Wpola 



fC2C^ 



(70) 



We examine equations (69) and ( |70[ ) to discover which one is physically 
significant. In the Schwarzschild limit {S -^ 0) we find that they coincide. But 
let us consider the weak-field limit (/ — > 00). Equation (69 1, Qpoiar = lb2bi , 
can be expanded in powers of /. And the terms with powers of l ^^ a nd lower 
approach zero as Z — > 00 to yield the asymptotic limit of equation ( 69 1 : 

(71) 



Cpolar 



e\ 



In a similar treatment of equation ( 70 ) one obtains: 

lc2c{^ = f+3 + e2-45^ 

which we can disregard as unphysical because it incorrectly implies that the 
spin of the KBH influences the trajectory of a test-particle at infinity. This 



(72) 



situation differs from that described by equation ( 54 1 , which does not describe 



a trajectory. Test c alcu lations confirm that equation (69 1 is the appropriate 



choice and equation (^70| can be disregarded since the values of Ilso obtained 



by solving equation (70 1 are less than (or equal to) the results from equation 



(69) {ceteris paribus). Equation (69) applies to all bound orbits, hence Qpoiar 
can be evaluated over a range I > lLSO.prograde\ but it applies only to orbits 
governed by Xl. In the SBH case (5 = 0), Bi = P - (3 + e^) f* and B2 = l^. 
Therefore 



Qpoiar — tit 6 O 

which has an asymptotic behaviour similar to that of Qx ■ 



(73) 
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4.5 Some Characteristics of the Carter Constant Formu- 
lae 

4.5.1 The last stable orbit on the abutment. 

In Figures ^ and dGJ) one observes that the functions for Qx and Qlso intersect 
at a single tangential point, which represents the value of I of an LSO that lies 
at the abutment described by X^. . The equation Qx — Qlso — (see equations 



(63) and (62 1 ) can be solved to yield: 

^(4 + Xi+Xo^i7), (74) 



''LSO ^ abutment -. (-, 



where 



^16 



^17 



(2I6 (1 + e) (5 + ef (1 - ef S^ + X1X2X3) + l2y/3Zr7, 

(1 + e) (5 + ef 

{l08(l + e)(5 + e)'(l~e)^54 

+X1X2X3 {^ - ef S^ 
-il + e){9-e')\l}, 

6e^ - Oe"* - 60e3 + llle^ + 342e + 441, 



6 



and 



Xo = e 



Xi = e^ + 3e2 + 3e + 33, 



X2 = e^ + 3e2 - 9e - 3, 



X3 = e^ + 3e^ - 21e - 39, 



X4 = e^ + 3e^ - 5e + 9. 



It is at this tangential point that the Qlso curve is split into two segments: 
the minus segment (I < I lso. abutment) that defines the values of Ilso (and 
Qlso) associated with inclined LSOs that are governed by X'^; and the plus 
segment (l > Ilso, abutment), which corresponds to inclined LSOs governed by 
X^. Consequently, the points beneath Qlso define only orbits governed by X^. 
Further, orbits with I and Q values that lie to the left of the minus segment of 
Qlso are undefined. Therefore the Qlso curve for I < Ilso, abutment and the 
Qx curve for I > Ilso. abutment define a curve along which the upper limit of Q 
is specified. 



Calculations of the value of Ilso. abutment (equation (74l) were performed 



for various eccentricities and two values of KBH spin (S = 0.5 and S = 0.99); 
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Table 1: Parameters calculated for the tangential intersection of the Qx and 
Qlso curves for an elliptical orbit about a KBH of spin S — 0.5 and S — 0.99. 
The I of the LSO at the maximum point of the Qlso curve is included for 
comparison. 









S = 0.50 








S = 0.99 




e 


Qlso 


i (deg.) 


''LSO. abutment 


lLSO,max 


Qlso 


i (deg.) 


'■LSO^ abutment 


lLSO,max 


0.0 


12.054503 


97.6 


6.067459 


6.067459 


12.203171 


104.4 


6.245480 


6.245480 


0.1 


12.102329 


97.2 


6.256148 


6.268429 


12.245340 


103.8 


6.407328 


6.449754 


0.2 


12.237084 


96.9 


6.445623 


6.467077 


12.366092 


103.3 


6.570627 


6.646723 


0.3 


12.449596 


96.7 


6.635958 


6.663706 


12.559674 


103.0 


6.735976 


6.836993 


0.4 


12.734670 


96.6 


6.827228 


6.858566 


12.823381 


102.8 


6.903945 


7.020995 


0.5 


13.090081 


96.5 


7.019514 


7.051863 


13.156912 


102.7 


7.075097 


7.199015 


0.6 


13.515997 


96.5 


7.212909 


7.243768 


13.562022 


102.7 


7.250004 


7.371204 


0.7 


14.014676 


96.6 


7.407517 


7.434425 


14.042359 


102.9 


7.429271 


7.537586 


0.8 


14.590357 


96.7 


7.603465 


7.623953 


14.603445 


103.1 


7.613543 


7.698050 


0.9 


15.249308 


96.9 


7.800900 


7.812451 


15.252772 


103.5 


7.803527 


7.852335 


1.0 


16.000000 


97.1 


8.000000 


8.000000 


16.000000 


103.9 


8.000000 


8.000000 



they are listed in table [ij We have also numerically calculated the values of / 
for the maximal point of the Qlso curve and listed them for comparison. For 
both circular and parabolic orbits, 



LSO, abutment 



= hso, 



irrespective of the value of S. Otherwise, 

'' L S O , abutment ^ ^LSO^max^ 



(75) 



(76) 



from which one may infer the curve that specifies the upper limit of Q is mono- 
tonically increasing with respect to / (with a point of inflection at Ilso, abutment 
for circular and parabolic orbits). The values of Ilso. abutment show that the 
upper limit of Q also increases monotonically with respect to e. Further, 
Ilso, abutment = 8.0 for a parabolic orbit, which equals the Ilso of a parabolic 
orbit about an SBH. 



4.5.2 The last stable polar orbit. 

The polar curve applies to polar orbits, which are governed by X"^; therefore, 
only the intersection of the Qpoiar curve with the minus segment of the Qlso 
curve (l < Ilso. abutment) needs to be considered. It is at this point that the 
^LSO of a polar orbit of arbitrary e is defined. It was found from numerical 
calculations of Ipoiar. lso (where the Qpoiar curve intersects the minus segment 
of the Qlso curve) and I value at the minimal point of Qpoiar that in the case 
of circular and parabolic orbits. 



polar, min 



polar, LSO 1 



(77) 
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otherwise, 

''polar, vain ^ ''polar, LSO- \'^} 

One may infer from this result that the value of Qpoiar is monotonically increas- 
ing with respect to I (and has a point of inflection at Ipoiar. lso for circular and 
parabolic orbits). 

4.5.3 The Schwarzschild limiting case. 

The analysis of these formulae in the case where S' ^- is an important test. An 



examination of the analytical fo rmu lae for Q switch (equation (51 1), Qx (equa- 



tion (63l), and Qpoiar (equation (69)), show that when S — 0, all three formulae 
equal 

Q = Pff-e^-s) . (79) 

In the Schwarzschild limit, we find that the abutment and the set of polar orbits 
approach one another, as required by the spherical symmetry of Schwarzschild 
spacetime. 



5 The Analysis of the Carter Constant for an 
Evolving Orbit 

5.1 Introduction 

We shall now perform an analysis of the evolution of Q in Kerr spacetime in the 
domain, which is defined by the three Q curves we have derived (QlsOi Qxj 



polar) in equations (62 1, (63 1, and (69 1. The behaviour of Q switch (equation 



(50)) will not be considered here, although it is important in guiding the choice 

of sign in ±a/X^. 

The three equations for Qx, Qlso, ^^nd Qpoiar define a map (see Figure p| 
from which one might infer the characteristics of a path followed by an inclined 
orbit as it evolves. These paths {Qpath) fall into two families: one governed 
by X^ and the other by X^. We conjecture that paths in the same family 
never cross; therefore, if Qpath reaches Qx, then it can do so only once. Let 
us concentrate on the behaviour of an evolving orbit at the abutment, which 
is where, Qpath, niay change from one family to the other. Therefore there are 
two modes to consider: 



Xl ^ Xl (80) 

and 

xl => xl. (81) 
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The mode represented by equation ( 80 1 corresponds to a rapid listing rate, where 
prograde orbits can cross Qpoiar and intersect Qx- The mode represented by 



equation (811 corresponds to orbits that: hst at a slow rate, have constant l, 



or exhibit decreasing t over time. And a prograde orbit cannot reach Qx since 
Qx lies on the retrograde side of Qpoiar- 

For this paper we will consider the evolution of a circular orbit (e = 0) 
because we wish to limit our initial analysis to the relationship between Q and 
/. Elliptical orbits will be treated in a forthcoming paper [38] . 



5.2 The Evolutionary Path in the Q - I Plane 

In Figure rol one may imagine a path, Qpath-, that starts at a large value of I as 
both Q and I monotonically decrease with respect to time. If the curve reaches 
Qx then it must intersect it tangentially (as the zone above the Qx curve is 

r2 



inaccessible). It is at that point that the orbit ceases to be governed by X^ and 



is governed by X^.. 

At the abutment, dQ/dl and dQ/de can be calculated analytically (see Ap- 
pendix B) irrespective of the model used to determine the radiation back re- 
action. Given an / = dl/dt and e = de/dt that have been derived according 
to some independent model, then according to a linear approximation one may 
define 



Q 



dQ-rdQ. 
dl de 



For a circular orbit (e — 0) it has been proven that e 
second term in equation (|82|) is zero. 



(82) 
|T5]; therefore, the 



5.2.1 A Preliminary Test at the Abutment. 

Because we can perform an a naly tical calculation of dQ/dl at the abutment, we 

can estimate I viz. equation ( 82 1 , 

-1 

(83) 



m 
Wt 



\dl 



if the 2PN Q flux (see equation (A. 3) in [55] (after equation (56) in [S])) is 
known, i.e. 



dQ 
dt 



2PN 






39 



(e) - [-I511 (e) + (^512 (e) - cos {i)Sg\^ (e)) P^/^ 



.9i3(e)-^'(5i4(e)-fsin2(i)) \l 



where the functions gg (e), g\f^ (e), gn (e), gi2 (e), 513 (e), and gu (e) are listed 
in Appendix B (the Carter constant, Q, has been normalised by dividing by 



(84) 
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Table 2: An estimate of (^M^/m) dl/dt, based on equation (83), for circular 
orbits. The values in parenthesis (which were originally calculated at t = tt/S) 
are taken from Hughes |1] and have been adjusted by dividing by sin (t). 

1 = 7 1 = 100 



5 = 0.05 -1.1827x10^1 
S = 0.95 -1.3747 x 10"^ 



(-1.2638 X 10^1) -1.2683 x 10"^ 
(-5.2540 X 10^2) -1.2679 x 10"'^ 



(-1.4637 X 10-'^) 
(-1.4553 X lO"'^) 



(mM) ). We performed test calculations on circular orbits of I = {7.0, 100.0} 
with KBH spin S = {0.05, 0.95}, which correspond to those used by Hughes 
|4] . In table [2] we compare our results with those of Hughes to find that they 
are reasonably consistent, with some deviation for S = 0.95 and / = 7.0. We 
have adjusted the results in Hughes by dividing them by sin (t) (where t = 
7r/3) so that they will correspond to our near-polar orbits; but this is only an 
approximation since t appears in other terms in equation ( 84 ) . 



5.3 The analysis of l on the abutment 

The analysis of u and its derivatives with respect to I and e will be treated 
in detail in a forthcoming paper |38| ; but it is appropriate to present a short 
preliminary exploration here. 

Unlike Q2PN (which appears in equation ( |84[ )), Qx contains no explicit 
variable, t; therefore, for circular orbits this property greatly simplifies the total 
derivative of Qx with respect to I since dQx/di = 0, i.e. 



(85) 



dQx 

dl 


dQx dQxdi 
dl di dl 




dQx 



dl 



Therefore equation ( 85 1 demonstrates that du/ dl is n ot constrained on the abut- 
ment in the same way as dQ/dl (see section 5.2.1). Consider equations (63 1 



(48), (44), and (42). They form a calculation sequence, which on the abutment 
creates a one to one mapping, Qx -^ t; otherwise, there are two possible values 
of 6 for a given value of Q. Thus dc/dl can be found either by numerical meth- 
ods or analytically [38 . In the remainder of this section we shall investigate the 
behaviour of dt/dl for orbits on the abutment. 



5.3.1 Numerical Analysis of l on the abutment 

We can numerically estimate the change of l with respect to / at the abutment 
by first finding the change in Q for an extrapolation of the evolving orbit's path 
(Qpath)- Because both Qpath = Qx and dQp^ff^/dl = dQx/dl at the point 
where Qpath intersects Qxi the equations of the second-order extrapolation of 
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Qpath at the abutment can be written as 

Qpath{l-Sl) = Qx-Sl 
Qpath (j+Spj = Qx+6l 



~dQx , SP d'^ Qpath 

2 dP 
~dQx , 5P d^ Qpath 



dl 



dP 



(86) 
(87) 



These equations are used to calculate t_ (associated with I ~ SI) and i+ {I + 61), 
where 61 is the small amount (10^'^^) by which we extrapolate from the value of 
I at which Qx and Qpath intersect. Equations (86) and (87) include the second 
derivative of Qpath with respect to I, which warrants further analysis. 

Because 61 is so small we used MATLAB, set to a precision of 256 digits, 
to perform the numerical analysis. From these extrapolated values of Qpath we 
obtain 



(k 
dl 



{l+-l^){261 



(88) 



N.B.: Since t increases as / decreases (with respect to time), di/dl ^ 0. 



5.3.2 The first-order extrapolation. 

Consider a first-order linear approximation, in which we drop the second deriva- 



tives in equations (86) and (87). The derivation of the corresponding change 



in L requires a sequence of calculations to be performed, which we will briefly 
outline. 

1. Specify the spin (S) of the KBH. 

2. Select the values of I and e for the point of intersection with the abutment. 
For this work, e = 0. 



3. Calculate Qx using equation ( 63 ) and dQx/dl using equation ( 105 1 (given 
in Appendix B). 



4. Calculate I lso, abutment usiug equation (74). It must be smaller than the 
value of / specified in point (p| otherwise the test-particle would be placed 
beyond the LSO. 

5. Calculate the values of Qpath {l~5l\ and Qpath \l + 6l\ according to a 
prescribed estimate or extrapolation at I. 



and we must use —\X^ if 



Cpath 



6. Calculate X^ and X+ (equation (48 1) from Qpath {l + 6l\ (We use —Jx\ 



I ±51) > Q switch)- 



7. Calculate X+ and X_ from Qpath (l — 61]. 
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Using equation ( 44 1 , calculate the orbital energies (E) for each of Qpath [l ~ SI 
Qpath {l + Sl 



and 



9. Now that the values of X+, X^,E. I are known; the value of l (t_ and i+) 
can be calculated viz. equation (42 1. We use the expression, L^ = X + SE. 



10. We can estimate di/dl viz. equation (88) 



We have p erfor med this sequence of first-order calculations for each mode (equa- 
tions (80 1 and ([8l| ) over a range of S and I values; and a representative set is 
shown in table I3j The slow and fast modes yield di/dl values that are of op- 
posite sign; and these results might suggest that the slow mode corresponds to 
orbits for which di/dt < 0. But let us first assess the validity of the first-order 
approximation by testing it in the Schwarzschild limit. The results for small S 
(table [3| demonstrate that this approximation is incomplete since it produces a 
non-zero result for 5^0, which is unphysical. On each side of the abutment, 
the value of Qpath is underestimated. This observation warrants the study of a 
more complete model that includes the higher derivatives of Qpath- Indeed, we 
have found that using dQx/dl alone cannot offer a sufficiently accurate math- 
ematical description of the orbital evolution at the abutment and warrants the 
development of second and higher derivatives of Qpath ■ One reasonably expects 
this numerical method to produce an accurate estimate of di/dl; but the transi- 
tion an orbit makes at the abutment from X'^ to Xj_ makes its behaviour more 
complicated. 



5.3.3 The second-order extrapolation. 



Equations ( 86 1 and ( 87 1 provide a second-order approximation of Qpath in the 



vicinity of the point of tangential intersection between Qx and Qpath- At this 
formative stag e of our work with the abutment, we will use d^Qx/dP (see 
equation (1071 in Appendix B) in place of d^Qpath/dl^ in equations (86) and 
(87) as an approximation. We repeated the ten-point sequence of calculations 



for the slow mode (equation (81 1), as outlined in section 5.3.2 The results of 
this numerical analysis are included in table [s] We observe that as ^ — )■ 0, 
di/dl — >■ 0, as required. Further, the second-order di/dl results for the slow 
mode represent listing orbits (i.e. di/dl < 0); therefore, the maximum list rate 
associated with paths that change from X'^ => X^ (slow) at the abutment will 
have an upper limit that corresponds to the minimal value of di/dl for that 
mode. 



5.3.4 The calculation of (di/dl 



Let us consider the use of the Qx curve itself to estimate the value of ( di/dl 
of an evolving orbit as it intersects the abutment. In this case, 

Qpath (i- Si) = Qx(l-Sl) (89) 
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Qpath{l + Spj = Qx[l + Spj, (90) 

where we have assumed that the path followed by the evolving orbit locally 
matches the Qx curve (equations ( 89 1 and ( 90 1 that are used in point Isl . This 

analysis yields the minimum value of ( db/dl J at I for a KBH spin. S, as 

0\ / min 

. For the slow mode, the rate of change of t can be no 

smaller. If Qpath deviates from Qx in its second and higher derivatives, then 
the actual value of db/dl will be greater than ( di/dl J 

V / min 

5.3.5 Analysis. 

In table 3 the values of ( di/dl J are in good agreement with the second-order 

I I V / min 

calculations; although a difference is evident for S — 10^^. 

We calculated (dL/dl) for various KBH spins (10^^^ ^ S ^ 0.99) over 

V / min 

a wide range of orbit sizes (10^ ^ I ^ 10^^). It was noted that the results for 
very large orbits were described well by an equation of the form 

'di 



^iSl' (91) 

ol. 



and that ki and p can be found by performing a least squares fit on I di/dl 
For orbits closer to the LSO, we find that 

^) = ~(^K,S + K,S')l'. (92) 

In Figure 7 ( dc/dl ) data for the range 10^ ^ I ^ 10^^ are shown on 

I I I V / min 

a log-log plot. By linear regression analysis, its asymptotic behaviour (/2.5) 
can be found (see table UJ) . In successive steps each power of I in the series 
can be derived as the higher powers are subtracted from the original numerical 
data-set. A linear relationship between IdL/dlj and S is found for /2.5 and 

V / min ^ 

/as; but for /4.5 and /as, which cover ranges of I closer to the LSO, an S^ 
term appears. The correlation coefficients (r^) of these regression analyses were 
better than 99.9999%. 

One may find the specific mode that applies at the abutment by comparing 
the results in the first line of table |4] with the weak-field radiation-reaction 
post-Newtonian results available in the literature. Consider the quotient of the 
formulae presented in equation (3.9) of Hughes [3] where t ^ 7r/2. 

itit^nh di 61 ;^~ 5 , ^ 

^^■^^^ = -. = --TT^Sr^. (93 

rweak a; 48 

An identical first-order result can also be derived from equation (4.3) in Ganz 



[7]. Equation (93) is of the same form as equation (91 1. Because —61/48 > —4.5 
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Table 4: The coefficients and powers of the series that describes ( di/dl ) 





Range 


Q 


Kl 


\ / rain 


h.b 


10^ < ^ < 1012 


-2.500 ±3 X lO"'^ 


4.500 ±5 X 10-13 




fs.b 


10^<[<10i2 


-3.500±7x 10"^ 


7.500±lx 10-^ 




/4.5 


10^<f<10i2 


-4.500 ±2 X lO"'^ 


31.5±4x 10-^ 


8.75±lx 10-'' 


h.b 


102<f<10i2 


-5.500 ±2 X lO"'^ 


122.70 ±2 X 10-^ 


-35.98±9x lO-'' 



Table 5: An estimate of (^M'^/m) di/dt, based on equation (94 1, for circular 
orbits. The values in parenthesis (which were originally calculated at t = 7r/3) 
are taken from Hughes 0] and have been adjusted by dividing by sin (t). 

f = 7 1 = 100 



S = 0.05 -9.4541 X 10"^ (-1.2557 x 10"^) -2.9301 x 10" 
S = 0.95 -1.8144 X 10-3 (-1.3941 x 10-^) -5.5681 x 10" 



" (-7.7291 X 10-12) 
1° (-1.3903 X 10-1°) 



in the weak-field radiation-reaction regime, one may consider X'^ => X^ to be 
the pertinent mode. Therefore (di/dl) describes the lower limit of di/dl 

mm 

for all I > I Lso, abutment- By numerical analysis, we found 

-_] = - (l22.7^ - 36^3) [-11/2 - (63/2^ + 35/4^3) f-^/^ 

-15/25f-^/2 - 9/2St^/^. (94) 

In Table Is] we compare the results of equation ( [94| with those of Hughes 
U; and although they differ, it is confirmed that the listing of an inclined 
orbit in a KBH system proceeds by the slow method. As before (Table l2| 
we adjust the results in |4] by dividing them by sin {tt/3) so that they will 
correspond approximately to our near-polar orbits. Although we use equation 



(42) to calculate the value of t in this work, and this differs from the formula 
used in pil |5] ; we recognise that they are sufficiently similar for us to make a 
general inference about the relative sizes of dt/dt in [¥ and those calculated 
here. 

Figure plshows the contours of constant Q on an / — i plane for circular orbits 
(e = 0) about a KBH of spin S = 0.99. One of the important features of dt/dl 
on the abutment is the suggestion of a coordinate singularity (di/dl — > oo); but 
this is for the specialized case in which the orbit evolves with a constant value o f 
Q. It has been confirmed that dQ/dl > on the abutment (see Section 5.2.11; 



hence, such a singularity for the di/dl of an evolving orbit is not physically 
manifested. The arrows labelled (a), (b), (c), and (d) show some important 
examples of how di/dl can vary at the abutment. One observes that di/dl 
is not uniquely fixed by dQ/dl. Nevertheless, Figure ^ provides an important 
picture of the behaviour of t as the orbit tangentially intersects the abutment; 
and it warrants further study. 
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6 Conclusions 

In our study of inclined elliptical orbits about a Kerr black hole (KBH), we 
found that the minus form of X^ (X — L^ — SE) is shifted away from the polar 
orbit position to encompass near-polar retrograde orbits. The abutment (which 
is a set of orbits that lie at the junction between the minus and plus forms of 
Xj.) is shifted the greatest near to the LSO of the KBH and asymptotically 
becomes more polar with increasing latus rectum (I) . 

We developed a set of analytical formulae that characterise the behaviour 
of the Carter constant (Q) at the last stable orbit (LSO), abutment, and polar 
orbit. Further, the curves that describe Q for an LSO and Q at the abutment 
(between the minus and plus forms of X|.) intersect at a single tangential point, 
for which we derived an analytical formula. From these equations one can 
define the domain of Q for an evolving orbit (Qpath)- The two families of curves 
defined by Qpath are governed by either X"^ or X^ ; and the curves within each 
family never cross. Therefore, at the abutment, Qpath can either change from 
X'^ ^ Xl_ or from X'^ ^ X£ . This result aids in the investigation of the listing 
of an orbit at the abutment. 

We have used the abutment as an analytical and numerical laboratory for 
the study of the evolution of Q for inclined circular orbits. The first derivative 
of Qx with respect to / (dQx/dl) allows us to test the consistency of 2PN Q 
fiuxes with estimated values of dl/dt. Further, by converting Q to the angle of 
orbital inclination (t), it was possible to calculate the minimum rate of change 

of t with respect to I, (di/dl) , independently of a radiation back reaction 

V / min 

model. Comparison with published weak-field post-Newtonian results show that 
the X'^ => X^ mode applies; and this mode must apply to the entire range of 

orbit size, I ^ I L so, abutment- 

Although Qx and dQx/dl are important, the higher derivatives also display 
critical behaviour. The second derivatives of Qpath warrant more study as it will 
improve our understanding to the effect of the radiation back reaction on the 
listing behaviour of the orbit. The analysis of elliptical orbits at the abutment 
will introduce new elements to the listing behaviour, which arise from the first 
derivative of Qx and the second and higher derivatives of Qpath, both with 
respect to e. Such a result might be valuable in our understanding of current 
and future radiation back reaction models. 
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Appendix A: Terms 



Table 6: Orbital Parameters 

Parameter Symbol 



Normalised Symbol* 



Test-Particle Mass 

Mass of MBH (typically 10^ M©) 

Orbital Radius 

Semi-Major Axis 

Latus Rectum 

Proper Time 

Orbital Energy 

Effective Potential Energy 

Spin Angular Momentum (KBH) 

Orbital Angular Momentum (z component) 

Orbital Angular Momentum (9 component) 



Governs prograde, polar, and retrograde orbits up to the abutment 
Governs the retrograde orbits beyond the abutment 

*We set the speed of light and gravitational constant to unity (i.e. c = 1 and G = 1); 
therefore, mass-energy is in units of time (seconds) . 



m 


- 


M 


- 


r 


R = rM~^ 


a 


A ^ aM-^ 


I 


I = IM-^ 


T 


f = tM-^ 


E 


E = Em-^ 


V 


V = Vm-^ 


J 


S=\3\M-^ 


iz 


L, = L,{mM)~^ 


Le 


Lg = LeimM)-' 


S 


t = EM-2 


= i?2 + S-^cos^ (6) 


A 


A = AM-2 


le abutment 


J^ HL.-SE) 



Appendix B: Ancillary Equations 
The Kerr metric and its Inverse 

The inverse Kerr metric expressed in the Boyer-Lindquist coordinate system: 

-2M 



0/9^ 

M^(it:^+5^)^-M^S^A 



9al3 



Kerr 



A-M^5^sin^(9) 







P" 



^ 
A " 



-2M 



M''SRsm''{8) 







'(e) 



To simplify the presentation of the metric, we define the parameter: 



sin^ (e) 



.(95) 



M^ R^ + 5^cos^ [6] 



(96) 
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The inverse Kerr metric is: 



S"f 



(97) 



(E)-i 











-2SR 



i?2-2i? + ^^0 

10 

-2Sfl g g -R^-2ii+S2cos2(9) 



m(r^-2R+S'^) 



The determinant of the Kerr metric was calculated to be. Det — — S^sin^ (9). 
Effective Potentials 



The effective potentials [TJ [Ml EHl IS] that appear in equations ([18]), ([19]), (|20l 
and (pTl): 



Vr (R) = T2 - a 



. _\ 2 



R^+{L,- SE 



Ve {e) = Q- cos^ [9) 



S'{l-E^ 



Ll 



sin^ (6*) 



(98) 



(99) 



14 



sin^ (B) ] A 



(100) 



with 



Vt 



■(sEsin^ie) 



L: 



R^ + S'^]T 



A 



T = e(r? + 5^) - U 
The Carter constant (Q) is normalised, 

Q = 



(mM) 
cos^ {9) Ll 
sin^ (6*) 



2 (';3^ r 2 

sin^ (6*) 



(101) 



1 r cOS^ {9) L ^ ^^2 ^ ^^^2 (g) ^2 (^2 „ j^2^ 



ll + cos^ [9) S^[l-E 



(102) 



32 



Ninth Order Polynomial in / for calculating Ilso 



(0 



f-4{S + e)l^ 



- (-36 - 2 S^e^ + 6 S^ - 24e + 4S^e - 4e^] r 

+4S'2(e + l)(-e2 + 2g-7)[^ 

-S^ (e + 1) (-S^e^ + 5 S^e^ -3S^e-9S^ + 16 Qe^ + 8 Qe + 24 q) ~f 

+8Q52 (e + 1)^ Ue^ + S^e~2e-iS^ + 6\ l^ 
+8 QS'^ (e + if {-2e^-e^ + 2Q- l) P - 16 S'^Q^ (e^ - 2 e + 3) (e + 1)^ P 



+16 S'^Q^ (2 e^ - 4 e + 3) (e + 1)'' / - 16 S^Q^ (e - 1)^ (e + 1 



0(103) 



The First and Second Derivatives of Qx 

Given: 

6 = f(f-e2-3), 

^3 = 2[-e^-3. 



From equation ( 63 1 the equation for Qx is 



2P 



ix 



a-v«-^^ 



(104) 



We obtain the following first and second derivatives of Qx with respect to / and 
with respect to e: 



di " 6 r^^ 2 y^f^e^ 



90 



2Ui-VC?"e2(0 



9Q 



X 



AeP 



de 



6(1-^2) 



-r(i 



r(l-e2)^i_^. 



2Ui-\/«-^: 



(106) 
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a^Q 



X 



dp 



2/ 
6 



21-" 



1 



(2reie3 - e 



d^Qx 



+4 ^6 
4/'3 
+ (l-e^)^ 



4 (??-6) 

1 (2[ei6 - 6 
"2 y^?^6 



2) .^^6ei + (3 + e2)- 



3/2 



y^F^ 



2Ui-\/ei^- 



-l(l-e^f + 2e^ 



V 



(;2([-3(l + e2))+4f52(3e2-l 



+86^ -f(l-e2) + 



v/eF^6 



+2 (5e2 + 1) Ui - V^' - ^■■ 



The 2PN flux for Q 

According to equation (A. 3) in [5T (after equation (56) in [F): 



(108) 



dQ 
dt 



2PN 






5 M2 



ff9 



(e) - Pi^n (e) + (^512 (e) - cos (0^5?o (e)) ^^/^ 



.9i3(e)-^'(.9i4(e)-fsin2(i)) )? 



(109) 



where 



59 = 1 + o^ 



510 = 


61 91 2 461 

= 8 + 4^+64^ 


511 = 


1247 425 2 
336 ' 336 "" 
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and 



. 97 2 
512 = '^+y^ 



44711 302893 ^ 
^^^ ~ 9072 ^ 6048 ^ 



33 95 2 
^" = 16 + T6^- 



N.B.: the Carter constant has been normahsed by dividing by (mM) . 

Appendix C: An explicit treatment of 6.9/ dr in the 
effective potential 

For the sake of completeness, we shall demonstrate a method for calculating a 
polynomial that describes dO /dr. 
Recall that: 

A = M^(r^-2R + sA (110) 

and 

S = /?2 = M^(R^ + S^cos^{e)\. (Ill) 

Therefore one may consider a normalised form of these equations: 



A = ^ - {R"-'2R + S' 



and 



N.B.: 



S=-^ = (i?2 + ^2^os2 (0) 

A A 

- f. (112) 

In working with the quantity, dr/dr, in [29 it did not matter about the 
division of the radial distance by the black hole mass, M, since the proper time, 
T, would also have been normalised in the same way i.e., 

dr d{M{f,)) 
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Indeed, one ought to consider the normaHsation of the proper time with respect 
to black hole mass. Although one could escape difficulties when only considering 
dr /dr, it is mandatory that the normalised proper time be explicitly considered 
when evaluating the quantity, dO/dr. The polar angle, 9, is already dimension- 
less, and as such cannot be normalised. Therefore 



d0 


d9 


d^ 


d{M{fj)) 




1 d0 




Mdf' 



(114) 

We may rewrite the 4-momentum in terms of d9/df in addition to dR/df, 
where X ^ L, - SE. 



P-r 



p^fdR^ ^p^(d9\ ^,.,__. (^^3^ 



By evaluating P ■ P and substituting the known relation for dR/df one obtains, 

= - sin^ (9) (l - E^\ R^ + 2 sin^ {9) R^ 

-(x^ + S^ + 2 XES - cos^ {9) S^ (l - eA ) R^ 

+2 (x^ + S^ cos^ {9) + 2X cos^ (9) ES - cos^ (9) S^ (l - eA\ R 
-{{X^ + S^ + 2 XES) cos^ (9) - CDS'* (9) (l - E^) S^) S^. 
We know from equation dsl that 

^^(^) - -(l-^2^i?4 + 2i?3 (117) 

- (x^ + S^ + 2 SEX + q] R^ + 2 {x'^ + Q)R- QS^ 
and we can substitute this expression into equation ( |116[ ) to simplify it thus: 

fl'sinH9)Ui^^^'\ (118) 

= 0x R'^ + Ox R^ 
-((x^ + S^ + 2XES) cos^ (9) - cos'' (9) S^ (l - E^) - sin^ (9) Q\ R^ 

+2 ((x^ + S^ + 2XES) cos^ (61) - cos'' (9) S^ (l - eA - sin^ {9) Q\ R 

-((x^ + S^ + 2 XEs) cos^ (9) - cos^ {9) S^ (l - &\ - sin^ [9) Q\ S^, 
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which factors to 



t'.in^e)U(^§^') (119) 
A ( sm^ (e) Q- (x^ + S^ + 2 XEs) cos^ (9) + cos^ (0) S^ (l - & 



and simpHfies to 



'^^\^\ = ■ 2/ ^ (sm^(6>)g- [X' + S^ + 2 XES] cos' {6)+ cos'' {9) S'[\-E 



dfy sin2(6') 



-2 / £08^(6*) ^ ^2 

''\sin\e) 



= Q-^ii ■ 2] J. ] -S' [l-E^] cos^ {0) . (120) 

One can equate: 

mMLg = m^ (— ] (121) 



and thus obtain: 



Le = ^^, (122) 

ar 



which viz. equation (118 1 confirms the relationship for the Carter constant. 
Further, we note that Lj. = Q at the orbital nodes and L^ = at the zenith or 
nadir of the orbit. It is in these respects that the Carter constant possesses a 
physical meaning for a KBH. 
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Figure 1: The relationship between L^ and pericentre, i?, for an SBH. Various 
values of Q are depicted. The polar LSO and abutment are superimposed (d, e). 
In table [7] some values of Rlso for this SBH system are listed. 



Table 7: Numerical values of Rlso estimated from Figure [T] for a circular LSO 
around an SBH. There is no circular LSO for Q > 12. 



Q 


label 


Rlso 


(Pro 


grade) 


label 


Rlso 


(Retrograde) 


0.000000 


a 




6.000 




A 




6.000 


5.000000 


b 




6.000 




B 




6.000 


10.000000 


c 




6.000 




C 




6.000 


12.000000 


d,e 




6.000 




D, E 




6.000 


13.000000 


f 








F 






16.000000 


9 








G 
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Figure 2: The relationship between L^ and pericentre, _R, for a KBH with 
S = 0.50. Various values of Q are depicted. The polar LSO (d) and the LSO 
at the abutment (e) are distinct. In table [s] some values of Rlso for this KBH 
system are hsted. 
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Table 8: Numerical values of Rlso estimated from Figure |2] for a circular LSO 
around a KBH of spin S = 0.50. There is no circular LSO for Q > 12.0545. 



Q 


label 


Rlso 


(Governed by 


x-i) 


label 


Rlso 


(Governed by X'^) 


0.000000 


a 




4.233 




A 




7.555 


5.000000 


b 




4.709 




B 




7.227 


10.000000 


c 




5.366 




C 




5.366 


11.828365 


d 




5.842 




D 




6.287 


12.054503 


e 




6.067 




E 




6.068 


13.000000 


f 








F 






16.000000 


9 








G 
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Figure 3: The relationship between L^ and pericentre, _R, for a KBH with 
S — 0.99. Various values of Q are depicted. The separation of the polar LSO 
(d) and the LSO at the abutment (e) is increased with the higher value of S. 
In table [9] some values of Rlso for this KBH system are listed. 
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Table 9: Numerical values of Rlso estimated from Figure [3] for a circular orbit 
around a KBH of spin S = 0.99. There is no circular LSO for Q > 12.203171 



Q 


label 


Rlso 


(Governed by 


x-i) 


label 


Rlso 


(Governed by X'^) 


0.000000 


a 




1.455 




A 




8.972 


5.000000 


b 




3.074 




B 




8.403 


10.000000 


c 




4.730 




C 




7.501 


11.252920 


d 




5.280 




D 




7.091 


12.203171 


e 




6.245 




E 




6.245 


13.000000 


f 








F 






16.000000 


9 








G 
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Figure 4: A plot of X'^ with respect to Q for a circular orbit {I — 6.25) 
about a KBH of spin S — 0.5. The slope of X„ can be assumed to have 
no discontinuities; therefore, the point at Q switch ~ 11.26 indicates that if 

Q > Q switch then X_ = -xlx'l 
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Figure 6: The three Q formulae derived in Section ffl define a map. In zone (A), 
only prograde orbits are found. And in zone (B), both prograde and retrograde 
orbits are found. Above the Qpoiar curve (zones (a) and (b) ) only retrograde 
orbits can exist. But in zone (a), the orbits are governed by X^; while in zone 
(b). they are governed by either X^ or X^. The points along Qx and Qlso 
mark the values of u. In this case the orbit is circular (e = 0) and the KBH spin 
is ^ = 0.99. 
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Figure 7: The values of 



(dL/dl) 

V / 7 



plotted for various values of S for circular 



orbits, where S = 0.99 and 10^ < l< 10^^ 
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Figure 8: Contours of constant Q in the I — t plane for a circular (e = 0) orbit 
about a KBH with spin S = 0.99. Polar orbits are indicated by the short-dashed 
line on the l-axis. The long-dashed curve corresponds to the abutment. Four 
curves (solid lines) of constant Q = {12.25, 13.0, 14.0, 15.0} are shown over 
a range of orbital inclination angles (90° ^ t ^ 115°). The segment of each 
curve that lies below the abutment is governed by X^; above the abutment, the 
segment is governed by X^ . At the points of intersection between the abutment 
and the curves of constant Q, dt/dl = oo, which suggests a singularity. The four 
arrows represent four tangential intersections on the abutment: (a) corresponds 
to the case where t is constant; (b) corresponds to the evolution of the orbit along 
the abutment; (c) represents the fast mode, and (d) the slow mode. N.B.: the 
four cases cannot occur together; they are shown on a single plot for illustrative 
purposes. 
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